x=-5:0.1:5;
y=sin(x/2);
s=(2065*cos(1/2)-3780*sin(1/2)).*x.^3+(2280*sin(1/2)-1245*cos(1/2)).*x;
sub=y-s
plot(x,y,x,s,x,sub);